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1 Introduction 



Exact solutions of non-trivial problems are always of great importance. They give a hint 
about the structure of real problems and also provide a laboratory for testing approximate 
methods. The non-relativistic many-body Calogero-Sutherland models 0, Q together 
with their supersymmetric extensions [|, ^ provide one of the most valuable sources 
of exact solutions to one-dimensional many-body quantum mechanical systems. The goal 
of this article is to try to uncover a hidden reason for the solvability of these models in 
order to answer the question what is special about them and whether there are other 
exactly-solvable or quasi-exactly-solvable many-body problems. 

The remarkable discovery of the solvability of the bosonic A^-body Calogero 0] and 
Sutherland |^ models was at the time a state-of-the-art achievement. Of course, this 
raised the question of whether there existed a regular procedure to generate these models. 
After a few years, it was found that the models are connected with the root systems 
of the Ai\i_i Lie algebras [0, and can be obtained from the Laplace-Beltrami operators 
defined on the symmetric spaces; this procedure was called 'the method of Hamiltonian 
reduction' [0, ^ ^ The Hamiltonian reduction method provides a regular basis for 

an explanation of the solvability of the Calogero and Sutherland models, at least for a 
selected set of values of the coupling constant (s). Considerations of the root systems of 
the other simple Lie algebras made it possible to find several additional families of exactly- 
solvable multi-dimensional Schrodinger equations and to prove the complete integrability 
of all these models. 

Recently, another explanation of the exact- solvability of bosonic many-body problems 
was presented, which was based on the finding that the eigenf unctions of the iV-body Atv-i 
Calogero-Sutherland models form a flag coinciding with the flag of the finite-dimensional 
representation spaces of the Lie algebra gl{N) |12[. It was shown that the Hamiltonians 
of the Calogero and Sutherland models are nothing but different non-linear elements of 
the universal enveloping algebra of a Borel subalgebra of the gl{N) algebra. Unlike the 
"method of reduction", in the second method the coupling constants appear as certain 
parameters fixing the element of the universal enveloping algebra. The second method 
can be used to explain the solvability of the above models for all allowed values of the 
coupling constant(s). In the present paper we show that the BCn, Bn, Cn and 
Calogero (rational) and Sutherland (trigonometric) Hamiltonians, for any values of the 
coupling constants, are second order polynomials in the generators of the gl{N+l) algebra. 
These results were made possible by exploiting the fact that the configuration space of 
the above quantum-mechanical systems can be parametrized by variables in which the 
permutation symmetry is already encoded. The most suitable variables are given by the 
elementary symmetric polynomials of the coordinates of the particles. In these variables 
the eigenfunctions have an especially simple form. 

Other approaches have also been used to study the Calogero and Sutherland models. 
Some of these have been directed towards obtaining closed expressions for the wave func- 



tions. The bosonic Calogero model was solved using an operator method in [13 



similar method was applied to the bosonic Sutherland model in |T^. The eigenfunc- 



tions of the BCn Sutherland model has been studied in ||T0|, |15[ (ground-state) and [|6 



'For a contemporary viewpoint on the Hamiltonian reduction method see, for example, 
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It is worth mentioning the remarkable observation that the wave functions of the bosonic 
Aj^ Sutherland models are in correspondence with singular vectors of ly-algebras (see, 
for example, |T7[|). 

Quite recently, supersymmetric extensions of the many-body A]y Calogero ^ and 
Sutherland models were constructed using the standard prescription of supersym- 
metrization. In the present article we derive the supersymmetric extensions of the BC^, 
Bn, Cn and Calogero and Sutherland models and show that their solvability can 
be explained by the existence of the hidden Lie superalgebra gl{N + 1\N). It is worth 
mentioning that a supersymmetric analogue of the Hamiltonian reduction method has 
not so far been constructed for these models. 

The paper is organized as follows: In the next Section we briefly review the bosonic 
many-body Aj\f Calogero and Sutherland models with emphasis on the property of exact 
solvability, and also set up our notation. In Section 3 we study the supersymmetric 
extension of the Calogero model and show that it is exactly solvable. In Section 4 the 
same analysis is carried out for the supersymmetric Sutherland model. Section 5 is devoted 
to the Calogero models connected with the Lie algebras BCn, B^, Cn and Djy as well 
as their supersymmetric extensions; it is shown that those models are exactly solvable. 
In Section 6 the BC^, Bjy, Cn and Dn Sutherland models and their supersymmetric 
extensions are explored and their exact solvability is established. The Conclusion contains 
a summary of the results obtained and a discussion of some possible directions for future 
investigations. In Appendix A we present a realization of the Lie algebra gl{N) and 
the Lie superalgebra gl{N\M) in terms of first order differential operators. Appendix B 
contains the Lie algebraic forms of the Hamiltonians discussed in this paper. Finally, 
in Appendix C we give some details on the derivation of the Lie algebraic forms of the 
Hamiltonians for the Calogero and Sutherland models. 



2 Bosonic (many-body) Ajv- 
land models (review) 



Calogero and Suther- 



In this section we briefly review, following |T2[, the exact-solvability of the bosonic many- 
body Calogero and Sutherland models or, more precisely, the An-i Calogero and Suther- 
land models^. 

The Calogero and Sutherland models describe systems of N identical particles situated 
on the line and the circle, respectively. The degrees of freedom of these models are 
parametrized by N real coordinates Xi. The Hamiltonian of the Calogero model is defined 
by 

92 , ,1 _ 1 



n 
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1 ^ 

2it 



dx'" 



I 2 2 



Kj 



(2.1) 



where g = z/(z/ — 1) > — ^ is the coupling constant and u is the harmonic oscillator 
frequency. For convenience only, we place the center-of-mass under the influence of the 



^For the sake of simplicity, in Sections 2-4 we refer to these models as the Calogero and Sutherland 
models. 
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harmonic oscillator force. The ground state eigenfunction is given by 

¥^\x) = A'^(x)e-"^ , (2.2) 

where A{x) = Yli<:j \xi — Xj\ is the Vandermonde determinant and = Y^i^h Q As has 
been shown by F. Calogero, any eigenstate of the system can be written in the form 

*(x) = ^i'\x)P,{x), (2.3) 

where Pc{x) is a certain polynomial in the Xj's, symmetric under the permutation of any 
two particles. We refer to these polynomials as Calogero polynomials. The operator having 
the Calogero polynomials as eigenfunctions is obtained from the Hamiltonian ( p.l| ) by a 
"gauge" rotation, by which we mean the similarity transformation 

hci = i^i'Y'ncM^ . (2.4) 
The Hamiltonian of the Sutherland model is defined by 

1 ^ g 1 

^t^idxl 4f^;sm'(i(xfc-x0) 

where g = z/(z/ — 1) > is the coupling constant^. The ground state eigenfunction is 
given by 

¥,'\x) = {A^'"^\x)r , (2.6) 

where A^*'"^\x) = Y\i<cj I sin(|(a;j — Xj)\ is a trigonometric analogue of the Vandermonde 
determinant. It was shown by B. Sutherland that similarly to the Calogero model any 



eigenfunction of the Hamiltonian (|2.5| ) can be written in the form 

vl/(x) = ¥^\x)Psien, (2.7) 
where P<j(e*^) is a certain polynomial in e*^-', symmetric under the permutation of any two 



particles. These polynomials are the so called Jack-Sutherland polynomials |18, |^ (for a 



general description, see, for example, [IT^ PD|). Performing a gauge rotation of T^suth with 



the gauge factor \E'q*^ we arrive at the operator 

/^suth = i^i'Y'Hsuti.'^i'^ , (2.8) 

which has the Jack- Sutherland polynomials as eigenfunctions. 

In order to study the internal dynamics we introduce the center-of-mass coordinate 
Y = Ylf=i Xj and the translation-invariant relative Perelomov coordinates 



Vi = 1 = 1,2,. ..,N , (2.9) 



^One should stress the point that to a fixed value of the coupling constant g there correspond two 
different values of the parameter v, namely, v = a and u — (1 — a), giving rise to two famili es o f 



eigenfunctions. Of course, a should be chosen in such a way as to mini mize the eigenvalue; then (2.2) 



corresponds to the ground state. The value v = {1 — a) inserted in ( p.2| ) describes the ground state (if 
normalizable) but of another family. If g = 0, this family comprises the states of negative parity with 
respect to permutations. 
^See footnote 7. 
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which obey the constraint J2i=i 2/i = 0. In order to incorporate the permutation symmetry 



and the translation invariance we consider two sets of coordinates |T2 

{Xi,X2 

and 



where 



(r,r„(x) =a„(y(x))| n = 2,3,...,iV), 

(xi, X2,... xn) ^ {e^^, Vn{x) = a„(e*^(^')) I n = 1, 2, . . . , (iV - 1)) 
akix) = 



Erf If ... 'T' 

il<i2<-<ik 



(2.10) 

(2.11) 
(2.12) 



are the elementary symmetric polynomials (see, for example, [^]). crk{y{x)) are thus 
the elementary symmetric polynomials with translation-invariant arguments]^. Then the 
following statement holds [p!2| : 



After separation of the center- of-mass, the operators /icai o,nd /isuth? when 
written in the coordinates r l \2.1Q ) and rj ( 2.11\ ), respectively, are quadratic 
polynomials in the generators of the Borel subalgebra of gl{N) spanned by the 
operators Jf^ and J~ realized as first order differential operators (see Appendix 
A, eq. 1^): 

N N N N 

h= ^ AijkiJijJli+ ^ BijkJijJk + ^ CijJij + ^DiJ~ , 

i,j,k,l=2 ij,k=2 i,j=2 i=2 

(see Appendix B). The operators J^j and J~ can be represented by triangular 
matrices and preserve the flag of spaces of inhomogeneous polynomials 



?14 



. . .V 



N 



where Vk = Tk and Vk = respectively, and k = 2,...N. The coupling 

constants g (see (2A ) and (2^)) appear only in the coefficients Cij and Di and 
are not related to the dimension of the representation of gl{N). Consequently, 
each Hamiltonian ^2. ^2.d[ ) has one or several infinite families of polynomial 
eigenf unctions. 

This statement leads to the important conclusion that the Calogero and Jack-Sutherland 
polynomials are intimately connected with the finite-dimensional irreducible representa- 
tions of the Lie algebra gl{N). It also provides a simple computational tool for deriving 
the explicit form of the Calogero and Jack-Sutherland polynomials. 



3 The super symmetric many-body Calogero model 

In order to proceed to the problem of the supersymmetric generalizations of (|2.1| ) and 
( p.5|) , let us first recall |Q that in a supersymmetric system of particles each bosonic 

^Due to the constraint X^i^i Hi — ^^'^ symmetric polynomials can be considered as defined in 
and then restricted to the hyperplane ^1^=1 Vi ~ ^- 
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degree of freedom Xi is accompanied by the fermionic variables ^9, and 'd\, which obey the 
standard anti-commutation rules {i^j, t^j} ={'(9j,'i9]} = and {-i^j, t?]} = Throughout 
the paper we use a concrete realization of this algebra: 

It is convenient to introduce the 'cmino' \1/ — the fermionic analogue of the center-of-mass 
coordinate as 

N 

^ = E^- (3-1) 

i=l 

which is the super-partner of Y. We can also define the fermionic analogue of the Perelo- 
mov coordinates ( |2.9| ) 

A. = - . (3.2) 

In order to construct the supersymmetric many-body Calogero model let us introduce the 
supercharges Q and Q"^ as defined in @]: 

where = Q^'^ = 0, Pk = ""^af^; the superpotential is 

W{xi,X2, . . . ,xn) = + ^^^og\xi - Xj\ . (3.4) 

^ i=l i<j 

Then the supersymmetric Hamiltonian TisCai = has the form 



1 



TV , g2 



2tt 



N 



d ^ V 



i=l i<cj i^i ^ji 



(3.5) 

where C = —}^vujN{N — 1) — \Nuj. The ground state eigenfunction remains the same (cf. 
( |2.2|) ) as for the bosonic many-body Calogero model ( p.l| ). It is easy to see that a gauge 
rotation of TisCai with the ground state wave function as a gauge factor affects only 
the bosonic part of the Hamiltonian. Defining 

Kci = -2(v|/(^y^7^3Cai^f^ (3.6) 

(cf. O), we obtain 



(3.7) 

This expression can be called the rationaHoim of the supersymmetric many-body Calogero 
model. Similarly to the bosonic case the eigenfunctions of the operator ( p.7| ) after sepa- 
ration of the center-of-mass remain polynomials but now in (y*, A*). These polynomials 
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are symmetric under the permutation of any pair of particles (x*,^*) ^ {x^,6^), and can 
be considered as the supersymmetric analogue of the Calogero polynomials. 

Now let us introduce permutation-symmetric variables. One can construct two sets of 
such variables: (i) a set analogous to the Newton polynomials 



N N 

Sn = Y.^1 , Pn = T.^^^r' (3.8) 
i=l i=l 

and (ii) a set analogous to the elementary symmetric polynomials 

(^k — 2^41^12 ■ ■ ■ -^ifc ; Cfc ~ 77 pTj (^i^Xi^ • " ' Xi^^ . (3.9) 

The variables s„ and p„ as well as 0"^ and are symmetric under the permutation, 
(a;\6'*) {x^,9^). However, the property aN+k = , (N+k = , = 1,2,... implies 
that the variables ak and (k, has the advantage that they avoid the difficulties associated 
with the overcompleteness of the basis, which plague the variables s„ and pn- Therefore, 
in the sequel we make use of the variables and (k- 

It is worth mentioning the following relations between the two sets of variables 

N CO (_-i\n+l 

Y.ak{x)t'^ = expE^^ s„(x)r], 

k=0 n=l 
N oo oo / -| \n+l 

Y.Ux)t' = (^(-ir+Vm(x,^)r)expE^^s„(x)r]. (3.10) 

k=0 m=l n=l 

A convenient way of succinctly writing these relations is obtained if one introduces the 
"superspace" coordinates 

Xi = cri + aCi, (j)i = Xi + adi, (3.11) 
where a is a Grassmann number, = 0; then 

N oo (T\n+1 

j:Xki<Pi)t' = eMT. -^^n{<P^)n , (3.12) 

A:=0 n=l ^ 



encodes both the relations ( |3.1CI| ). 

The "superspace" formulation makes it possible to write formulas in a more compact 
way. The "supercoordinate" (pi is clearly an = 1 superfield. However, one should 
stress that the models we study possess N = 2 supersymmetry, while the "superspace" 
formulation is not manifestly N = 2 supersymmetric. Nevertheless, later on we will see 
explicitly that the superspace (we drop the quotation marks hereafter) formulation turns 
out to be a useful computational aid. For example, one can rewrite the Hamiltonian hgCai 
in the superspace coordinates as follows 



TV r r 

/iscal = / dadaY^-^— - [ daW - [ daW . (3.13) 

J O0i OCPi J J 
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Here a is a Grassmann variable independent of a and 

^ — ^ — (3 14) 

50j dQi dxi ' 50j dOi dxi ' 

while 

where 0j = Xi + a9i. Notice that we have not introduced any superderivatives; ^ and ^ 
should be considered as "superfields" (since they do not contain derivatives with respect 
to t, a and/or a). 

Like we did in Section 2 we would like to separate the dynamics of the center-of-mass 
and the relative motion, and encode in relative coordinates the translation invariance 
and permutation symmetry of the system. In order to implement these requirements we 
introduce the translation-invariant and permutation-symmetric coordinates of the relative 
motion 

niy) = (^kivix)) , Kk{y,X) = Ckiyix),X{e)) , (3.16) 

for example, in the symmetric polynomials ( p.9|) the arguments (x's) are replaced by the 
translation-invariant y's ( |2.9| ) and the 6''s are replaced by the translation-invariant A's 
( p.2|) . Here Tq = 1, ri = and kq = ki = 0. The next step is to make the change of 
variables 

{xi, ei\t = l,2,...,N) {Y, Tk, ^,Kk\k = 2,3,...,N) . (3.17) 

in the Hamiltonian /isCai- In general, this is a tedious and cumbersome calculation. For- 
tunately, the superspace formulation allows us to avoid most of the tiresome algebraic 
calculations. Let us define 

1 ^ 

Then the expression for ipi = Ti + ani can be obtained from the relation 

A'" _ oo / -| \n 

J2 M4>i)t' = exp[^ ^-^Sn{4>^)n ■ (3.19) 

fc=0 n=l 

The derivative cj, , = ^ + satisfies 



6ij{a - p) = 6ij6{a - p) , (3.20) 



50j(a) 

where we have used that / da6{a — (3)f{a) = f{P). We now make the change of variables 

N 

(0,|2 = l,...iV)->(xi = E0fc; i;,\j = 2,...N) . 

k=l 

Under such a change of variables the derivatives transform asj^ 

(3.21) 



6(f)i{a) J 6(f)i{a) 5i)j{[3) 
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From now on we suppress the center-of-mass coordinate xi- 



7 



where 



(5i/'i(o) 



7^ + Ot^, and summation over j is assumed. To prove the relation (13.21 



one simply writes both sides of the equality in components. It is convenient to perform 
the change of variables in two steps: to write, 



5^i{a) 50^.(7) 5MP) 



(3.22) 



and then use the definition of 0j together with (|3.21|) which leads to ^^^^^j — i^ij ~ 
j^)6{a — 'j). Let us give as an example the expression for the Laplace operator in the new 
coordinates 

TV ^ ^ N N N 



1=1 



1 S 



(3.23) 



The above considerations imply that within the framework of the superspace formalism 
the calculations for the supersymmetric models follow closely the calculations carried out 
in the bosonic case |12[. In fact, for most of the terms the supersymmetric results can 
be obtained from the bosonic ones by judiciously replacing the bosonic coordinates by 
superspace coordinates. More details can be found in Appendix |C|. The derivation makes 
extensive use of the generating function ( |3.12|) . 

In the superspace coordinates the Laplace operator becomes (after reinserting the 
center-of-mass coordinate) 



TV 

/ dada 
k=i-' 



6 6 



dada 



dxi dxi i~^2 ^i^j 



TV e ^ X 



where 



and 



A 



i)i = Ti + aKi 
(iV-z + l)(j-l) 



5 _ d d 

dKi dTi 



(3.24) 
(3.25) 



N 



ipi-iipj-i + U - (3.26) 

Z>max(l,j— i) 



1>1 

1 N 



2N 



(3.27) 



i=2 



The final expression for the Laplace operator in components is given by 

TV ^2 N 



92 



92 



g2 Q2 Q2 

ATT ^ _|_ ATK. ^ _|_ AKT ^ _|_ Ah 



^ r (9 d 



i=2 



dKi 



(3.28) 
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where 



ATT 




(Ar_i + i)(j 


-1) 


— 


N 




ATK 






- 1) 




N 






— 


1>1 


- Ti 


AKT 


= 




- 11 




N 






+ 


l>l 


- Ti 


A KK 






-1) 




N 






+ 


2 ^ lKi+lK,j-2- 


-I ) 



i>max(l,j— i) 
Z>max(l J— i) 

Z>max(l,j— i) 

Z>max(l,j— j) 



Z>1 



{N -i + 2){N + 
N 

(N -i + 2){N + 
N 



-Ti-2 



(3.29) 



As expected, AJJ and BJ coincide with the expressions found in [jT2l for the bosonic 
Calogero modeL Let us note that the coefficient functions Aij and Bi are second and first 
order polynomials in Ti and Kj, respectively. 

Dropping the decoupled center-of-mass terms one can show that in superspace coor- 
dinates the remaining part of the Hamiltonian h^Jcl^ which describes the relative motion, 
becomes 



^sCai — J dada X 
where Aij is given in ( |3.26|) and 



6 S 

i,j=2 ^^j 



daW - / daW 



(3.30) 



TV r -, -, AT r 

Finally, in component form /igcii looks like 



(rel) 
sCal 



N 

E 



g2 Q2 Q2 Q2 



dKjdKi 



. 1 



N 



V 



i=2 

N r g 



d d 

T"j-2 7^ \- l^i~2Tr- 

OTi OKi 



i=2 



d 

. OTj OKj 



(3.32) 
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where AJT Aj-", A^T A^-^, BJ, are given in (|3.29|) . This form can be called the algeh 



raze 



form of the supersymmetric many-body Calogero Hamiltonian. 

It turns out that we are able to rewrite the supersymmetric Calogero Hamiltonian 
h^scJi in terms of generators of the Borel subalgebra of the algebra gl{N\N — 1) spanned 
by first order differential operators in the variables (r/j, Hk), k = 2, . . . N (see Appendix 
A. 2). The result is given in Appendix B, eq. ( [B.2|) . Since the expression (|B.2| ) contains no 
positive- root generators J^^, Q^, then in accordance with the general definition given in [|l| 
we conclude that the supersymmetric many-body Calogero model ( WM is exactly solvable. 
The existence of the representation ( p.2|) proves that there are infinitely-many eigenfunc- 
tions of the operator (3.7) having the form of polynomials in the variables {rk, Kk)- It also 
implies that eigenfunctions of the supersymmetric many-body Calogero model ( |3.5| ) have 
a factorizable form being the product of the ground-state eigenfunction multiplied by a 
polynomial in the variables {Tk,Kk)- These polynomials are related to finite-dimensional 
irreducible representations of the algebra gl{N\N — 1) in the realization (A. 3). 



4 The supersymmetric many-body Sutherland model 

In this section we analyze the exact-solvability of the supersymmetric Sutherland model. 
In order to define the supersymmetric extension of the bosonic many-body Sutherland 



model ( |2.5| ) we use the standard prescription of supersymmetrization already used in the 



previous Section. Let us take the supercharges Q,Q^ (see ( ^731) ) with a superpotential 

W{xi, X2, . . . , xn) = i^Elog I sin(^(x, - Xj))\ = i/log A^*"^) , (4.1) 

i<j 

(cf. ( ^^) ) and then construct the supersymmetric Hamiltonian TisCai = \{Q,Q^- After 
carrying out the calculation the Hamiltonian of the supersymmetric many-body Suther- 
land model emerges [§]: 



1^ 1 ^ z/ 



'^t^i^^l ' 4f^;Sin2(l(xfe-a;0) 



C , (4.2) 



where C = —N{N'^ — l)i^^/3. The ground state wave function remains the same as in the 
bosonic case, "^q^ = (A'^*"^^(x))'', where A(*"3)(a;) = llk<i I sin(i(a;fc — xi)\. Introducing 
the gauge-rotated Hamiltonian /isSuth = —2 {¥o^)-^n¥Q\ we get 



TV 



1 ,,,9 d ^ 1 (Ok-ei) , d d 



k=l ^-^k k<l 

(4.3) 

It is worth mentioning that if the operator ( [4.3|) is written in the coordinates e*^^, it 
appears in its rational form, which is a supersymmetric generalization of the rational 
form of the Sutherland Hamiltonian (see 0). 

One can show that like in the bosonic case (cf. (|2.8| )) after separation of the center-of- 
mass motion the eigenfunctions of the operator ( |4.3| ) remain polynomials but now in the 
coordinates (e*^''',Afc). These polynomials are symmetric under the permutation of any 
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pair of particles (xj,^,) {xj,9j) and can be considered as a supersymmetric analogue 
of the Jack-Sutherland polynomials. 

In the superspace coordinates (|3.11|) , (|3.14|) the Hamiltonian /isSuth can be written 

r _ ^ 6 6 r r _ - 

/issuth = J dada ^^j- + 1^ J daW + v j daW , (4.4) 



k=l 



where 

^i<P) = lT.^ot{li<P,-<P,)){^-^) , W = W(0). (4.5) 

k<l 



Next, we introduce the new variables (cf. (|3.16|) ) 

Tin + Oil^n = cr„(exp[i(?/fc + ckAfc)]) , n = l,...,N -1 , (4.6) 

and 

cn + aCiv = (TN{exp[i{xk + a6k)]) , n = N , (4.7) 

where = 0, is defined in ( ^l9|) and ?/j,Ai are given in ( p.9|) , ( ^^) , respectively. 
Furthermore, 

A'' oo /'_]^^n+l 

n 



5: r/.t'^ = exp[^ ^^—s^{e^y>^)r] 

k=0 n=l 

n 



J2^^kt' = (E(-l)"^Vm(e^'N V^'=)t™)exp[5: '-^—s^{e'y'')r] , (4. 

k=0 m=l n=l 

(cf. (|3.1CI| )), where we have set rj^ = 1 and fiN = 0. 

Following the same procedure as in Section 3 one can rewrite the Laplace operator in 
the new coordinates (^l6| ), (|4.71 ): 



N f)2 a N-1 

B1— + 5f — 



Q2 Q2 Q2 Q2 
— r /iii — t: h -r h /i,-,; 



+ E 



(4.9) 



where 










/IW — 

^ij - 


N 


-ViVj + E 


i - 2l)r]i+ir]j_i 






(>max(l,j— i) 






^ = 


j{N-t 
N 




- 2l)r]i+ij^j^i - 






i>max(l,j— i) 




i>i 


^ = 


N 




- 2l)fii+ir]j^i + 


Y.llfJ'i+i-iVj-i-i -Vi+i 




i>max(l,j— i) 




1>1 




j{N-t 
N 




- 2l)fii+ifij_i H 


- 2 ^ lfii+i_iHj_i+i , 




/>max(l,jr— i) 






= 


i{N-i) 
N 










i{N-i) 


-fJ'i ■ 






Bi' = 


N 







(4.10) 
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As expected, the coefficients A^^ and B]' coincide with the expressions found in [12] for the 



bosonic Sutherland model. Furthermore, the coefficient functions Aij and Bi are second 
and ffist order polynomials in rji and fii, respectively. Similarly to what happens in the 
bosonic case the Laplace operator ( [4.9| ) possesses infinitely-many polynomial eigenfunc- 
tions ||12|. These have the form of supersymmetric analogies of the Bethe-ansatz wave 



functions. 

After omitting the center-of-mass motion associated with and (n, the final expres- 
sion for the supersymmetric many-body Sutherland Hamiltonian /isSuth takes the form (cf. 
Section 3, eq.( |3.32|) ) 



h 



(rcl) 
sSuth 



iV-l 

E 



Am 



dr]idi]j drjidfij 



drijdfii 



N-l 



+ (t7 + ^)E*(^ 

i=i 



d d 
. or]i dfii. 



(4.11) 



where A^j , A^j" , Af^ , Afj' , B^ , B-" are given in (|4.10| ). This expression can be called the 
algebraic form of the supersymmetric many-body Sutherland Hamiltonian. 

In superspace coordinates the Hamiltonian /isSuth describing the relative motion can 
be written 



/■ A^— 1 6 6 f I" 
- /^SL = j dada E -4^.^^ + JdaW + J 



where 



i=l 



5 d 
6% dni drii 



d 



(cf.(ra)), while 
jiN- 



A 



N 



l>ma,x{l,j—i) 1>1 



(4.12) 

(4.13) 
(4.14) 

iji+i-iijjj-i+i] , 



(4.15) 

(cf. (|3.26| )). The method of the calculation of the coefficients Aij is similar to that 
presented in App.C. 

It turns out that the Hamiltonian governing the relative motion can be rewritten 
in terms of the generators of a Borel subalgebra of the Lie superalgebra gl{N\N — 1) 
(see Appendix A. 2) after substituting Kj — and Xj — r/j-i. The final result is 
given in Appendix B, eq. ( |B.4| ). Since the expression ( |B.4[ ) contains no positive-root 



generators Tj , then in accordance with the general definition ||1|] we conclude that 
the supersjTiimetric many-body Sutherland model (|4.2| ) is exactly solvable. The existence 
of the representation ( p.4| ) proves that there are infinitely-many eigenfunctions of the 
operator (4.3) having the form of polynomials in the variables (77^, fik)- It also leads to the 
conclusion that eigenfunctions of the supersymmetric many-body Sutherland model ( [4.2| ) 
have a factorizable form being the product of the ground-state eigenfunction multiplied by 
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a polynomial in the variables (77^, Hk)- These polynomials are related to finite-dimensional 
irreducible representations of the algebra gl{N\N — 1) in the realization (A. 3). 

To conclude, the supersymmetric many-body Sutherland models possess the algebraic 
form ( |4.11| ) and also the Lie-algebraic form ( |B.4| ) represented by second-order polynomials 



in the generators of the of the algebra gl{N\N — 1) with certain coefficients. 

5 The BCn, B^^ Cn and Calogero models and their 
supersymmetric extensions 

It is well-known that there is a deep connection between the integrable Calogero-Sutherland 
systems and Lie algebras (see discussion in the Introduction). This connection is explicitly 
realized in the Hamiltonian reduction method [0, |^ (see also |P, |lTl). In particular, 
the celebrated many-body Calogero and Sutherland models discussed in Section 2 are 
related to the Lie algebras v4jv_i and are therefore known as the Calogero and Suther- 
land systems of Ajv_i-type. A natural question emerges: what are the integrable systems 
corresponding to the other simple Lie algebras like Bj^f, Cn, Dn etc. appearing in the 
Hamiltonian reduction method. The answer and a complete classification of these systems 
is given in [0 (for a review, see |P, |l^). In the present paper we focus on the quantum 
Calogero (rational) and Sutherland (trigonometric) systems of the BCn, Bjy, Cn and 
types leaving for a future publication the case of the exceptional algebras. 



5.1 The bosonic case 



Unlike the Aiy_i many-body Calogero model, the BC^, B^, Cn and Dn Calogero models 
are not translation invariant, and describe systems with boundaries. However, they are 
still permutation invariant with respect to the interchange of any pair of coordinates. The 
configuration space of these models is {a;j|xj > 0; j < : < x^}. In this section we will 
show that like the An-i many-body Calogero system, all the BCn, Bn, Cn, Dn quantum 
Calogero systems are not only completely integrable, but also exactly-solvable possessing 
a hidden algebra gl{N + 1). 



It is well known that the Hamiltonians of the BCn, 
given by (see [0) 



Bn,Cn Calogero models are all 



1 ^ 

^ i=l i<j L 



N 



(5.1) 



where g = v{v — 1) and g2 = 1^2(1^2 — 1)- When the coupling constant g2 tends to zero, 

(c) 

the Hamiltonian 'Hbcd degenerates to the Hamiltonian of the Dn Calogero model. The 
ground state eigenfunction of the Hamiltonian (|5.1|) is given by 



N 



i<j 



1=1 



2 Z^i=l -^i 



(5.2) 



(cf.(il)). 
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One should stress the point|^ that to fixed values of the coupling constants g2 there 
corresponds two different values of the parameter i^(z/2): 1^(1^2) = a(a2) and (1 — a(a2)) 
giving rise to four families of eigenfunctions. Of course, a{a2) should be chosen in such a 
way as to minimize the eigenvalue and then ( p.2| ) corresponds to the ground state. The 
ground state can be denoted as (a, 02) • The other values z/(z/2) inserted in ( ^.2|) describes 
the ground states (if normalizable) of three ot/ier families: (1— a, 02), (a, 1— 02), (1— 1 — 
02) • If 92 = 0, the sectors (a, 1 — 0:2), (1 — a, 1 — correspond to totally antisymmetric 

states with respect to the refiections Xj < Xj. In what follows, for the sake of simplicity, 

we call ( |5.2| ) the ground state, of course keeping in mind the above discussion. 

Following the same approach as in Sections 2-4, we make a gauge rotation of (|5.1| ) 
with the ground-state eigenfunction as the gauge factor; h^scD = — 2\E'o ^7-^\E'o. A straight- 
forward calculation gives (we omit an additive constant) 



N 



1 



N N 



=1 ^« i=l 



We call this operator the rational form of the BCn, Bj^, Cm and Calogero systems. 

The Hamiltonian ( |5.1| ) as well as the operator ( |5.3| ) is invariant under the permutation 
of any pair of the coordinates Xi ^ Xj and also under the action of the refiection group 
Z®^: Xi —Xi, i = l,2,...,iV. An infinite set of eigenfunctions of (|5.3| ) are given 
as polynomials in the x's and are classified by the representations of the refiection group 
Ti®^ . For the sake of simplicity, we consider in what follows only the eigenfunctions of the 
Hamiltonian ( [5. 3D , which are totally symmetric under the refiection and permutation group 
actions. The totally (anti)symmetric- (anti)symmetric eigeniunctions under the refiection 
and/or permutation group actions are reproduced by a change of parameters u, 1/2 in ( [5.11 ). 
This fact implies that it is sufficient to look for eigenfunctions to the operator ( |5.3|) of the 
form 

^^'^^=F{xlxl...,xl) , (5.4) 
Now let us construct variables which encode the permutation invariance of the system. 



Analogously to what was done for the A^^i Calogero model |T2[ we use as new variables 
the elementary symmetric polynomials ak (see (2.12)) but with xf as arguments: 



N 



J2 &k{xf)t'' = exp 



k=0 



' °° ( — } 

E ^^S2n(x.)r 
.n=l n 



(5.5) 



(cf. (2.10), ( ^.lOD ). Finally, in the new variables the Hamiltonian h^BCD becomes 

h'^l-o=j:A,^^ + j:B,£ (5.6) 

where (cf.(3.28)) 

^ii = 4^(2/ + 1+ j -i)a,_,_iaj+i 

B, = 2[l + U2 + 2u{N - 3)]{N - j + l)a,_i - iujdj , (5.7) 

^^See also the discussion in footnote 7. 
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and a"fc = 0, when k > N or k < 0. The method of the calculation of the coefficients 
Aij,Bj is similar to that presented in App.C. This expression can be called the algebraic 
form of the BC^, B^,Cn and Calogero Hamiltonians. 

Similarly to what happened in all previously discussed cases the coefficients Aij,Bj 
are polynomials of second and first degree in a^, respectively. Hence, there exists a 
Lie algebraic form of the Hamiltonian and h^BCD can be written in terms of generators 
of a Borel subalgebra of gl{N + 1) (see Appendix A.l) as in Section 2. The result is 
presented in Appendix B, eq.(B.5). Then in accordance with the general definition given 
in we can conclude that all BCn,Bn, Cn,Dn Calogero models (|5.1|) are exactly 



solvable. The existence of the representation (|5.6| ) proves that there are infinitely-many 
eigenfunctions of the operator (5.3) having the form of polynomials in the variables dk- 
It also implies that totally (anti) symmetric- (anti) symmetric eigenfunctions with respect 
to permutations and refiections of the BCn, Bn,Cn and Calogero models (|5.1| ) have 
a factorizable form being the product of the ground-state eigenfunction (5.2) multiplied 
by a polynomial in the variables cr^. These polynomials are related to finite-dimensional 
irreducible representations of the Lie algebra gl{N + 1) in the realization (A.l) and can 
be called the BCn Calogero polynomials. 



5.2 The supersymmetric extensions 



The BCj\f, B]\r,CN and D^r Calogero models discussed above have also natural supersym- 
metric extensions. These can be constructed in a straightforward way. Let us consider 
the supercharges (|3.3D with the superpotential W given by 



N 



W = z/^(log \xi ~Xj\+ log \xi + Xj\) + z/2^1og |a;i| - ^ 



N 



(5. 



i<j 



i=l 



1=1 



(cf.(|3^)). After some calculations the supersymmetric Hamiltonian Ti-'f^cD = ^iQ^Q^} 
of the supersymmetric BC^, BN,C]y, Djy Calogero models emerges: 



l^sBCD 



1 ^ 



2 2 



^ i=l -^i 



N Q 



V 



1 



i2/ T tJj /J 



V- 1 



i<3 



■ (Xj -|- Xj) 



+ c 



(5.9) 



where the constant C = —2N[i'{N — 1) — z/2 — 1]. Again, like for the other supersymmetric 
extensions the ground-state eigenfunction remains the same as in the bosonic case and is 
given by (|5.4| ). Making the gauge rotation of the Hamiltonian (|5.9|) : 

^ScD = -2^0 ^7i:ScD^o, we get 



''sBCD 



N Q2 

i=l ^-^i i<j 



d 



d 



d 



Xi — Xj dxi dxj {xi — XjY dOi dOj' 
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+ 

i<j 
N 

- 2^E 

i=l 



1 , d d , (9i + 9i) , d d . 
( \ ) — — —{ \ 

Xi + Xj dxi dxj {xi + Xj^ 89 i 89. j' 



Xi 



d_ 

dxi 



_d_ 



N 



i=l 



1 d _ 9^d_ 

oc j doc y X n d0 j 



(5.10) 



We call this operator the rational form of the supersymmetric BCn, -Btv, Cn Calogero 
systems. 

The superspace coordinates Xi which are invariant under the permutation and parity 
transformations are defined through the relation 



TV 



T.Xk{(l)l)t'' = exp 

A;=0 



.n=l 



n 



(5.11) 



dAxi) + 



(cf.(5.5) and (3.19)), where (pi = Xi + a9i as in Sections 3 and 4, and Xfc(0f 
a(k{xj,'2xi9i). Here a and ( are the elementary symmetric polynomials defined in 
but of the new arguments: x ^ x'^,9 —>■ 2x9. Using the same technique as in Sections 3-4, 
one can find the representation of the supersymmetric models directly from the analogous 
representation for the bosonic cases (for details see App.C). Finally, in the superspace 
coordinates (|5.11|) the Hamiltonian h^^BCD ^^e form 



h. 



(c) 

sBCD 



i,j=l j=l 



d_ 



(5.12) 



where 



>Aij 



2^ {2(2/ + 1 + j - i)Xi-i-iXj+i - l[Xj+iXi-i-i - Xj+iXi-i-i 



l>0 



+ Xi+l-lXj-l — Xi+l-lXj-l\ 

Bj = 2[l + iy, + 2u{N-MN-j + l)xj-i-4u;jXj. 
For completeness we give the form of h^sBCD components 



(5.13) 



\BCD 



N r a2 - R2 - - Q2 



E 



daidQ 



N 



Y\2[l + U2 + 2uiN-Mj-N + l) 
d r d 



dCidaj 
d 



d 



(5.14) 



where 



Af/ = 4 5](2Z + 1+J- t)a,_i_,a,+i , 

l>0 

^ij = 2 E [2(2^ + 1 + j - ^)^^-l-lCj+l -^(^j-i-iCi+i - S^i+iCj-i-i + ^j-iCt+i-i - cri+i_iCj_/] 

Z>0 
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2 ^ [2(2/ + 1 + J - + l[aj^ 



l>0 



4 J2 [(2/ + 1+J- t)Q-i~iCHi - K^+lC: 



i+i^j-i-i 



Ki+l-lCj-l 



l>0 



(5.15) 



/^From the expression ( [5. 141 ) with the coefficients ( [5. 15] ) it follows that h^^BCD can be 
written as a quadratic polynomial in the generators of a Borel subalgebra of gl{N + 1|A^) 
and, hence, the -BCat Calogero model is exactly solvable. 

The existence of the representation of ( |5.14| ) in terms of the generators of a Borel 
subalgebra of gl{N +1\N) (see (B.6)), proves that the operator h^^BCD has infinitely-many 
eigenfunctions having the form of polynomials. These polynomials are related to finite- 
dimensional irreducible representations of the algebra gl{N +1\N) in the realization (A. 3) 
and can be called the supersymmetric BCn Calogero polynomials. 

So the supersymmetric BCn, Bn,Cn and Dn Calogero models possess the algebraic 
form and also the Lie-algebraic form (B.6) being represented by second-order polynomials 
in the generators of the of the algebra gl{N + 1|A^) with certain coefficients. 



6 The BCn^Bn, Cn and Sutherland models and 
their supersymmetric extensions 

Similarly to what was done in Section 5 for the BC^, B]\r,CN and Dj^ Calogero models 
the present section is devoted to a consideration of the BC^, B^, Cn and D]\f Sutherland 
models. 



6.1 The bosonic case 

The Hamiltonians for BCn, Bn,Cn, Dn Sutherland models are special cases of the general 
BCn Hamiltonian [f^l 



H 



BCD - 2 ^ 4 ^ 



1 



sin2(|(xj — Xj)) ~'~ sin2(|(xj + Xj)) 



92 \ - 

+ tE- 



N 



4 — sm-(x,; 



+-E- 

8 fesin 



N 



(6.i; 



where g = v{v — 1), g2 = 1^2(1^2 — 1) and g^ = 2/3(1/3 + 2z/2 — 1). /,From the general 
Hamiltonian the B^, Cn and Dn cases are obtained as follows: 



-Bat case: 1^2 = 0, 
Cn case: 1/3 = 0, and 
Dn case: 1^2 = ^3 = 0. 
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The ground state wave function is [|T^, 15 



n 

i<j 



sin(-(xi 



N 

n 



sm X,- 



\'^2 I 



sm 



(6.2) 



Again, one should emphasize 0that to a fixed value of the coupling constant g{g2)[g3] 
there corresponds two different values of the parameter 1/(1/2) [i^s] : '^('/'2)[^'3] = tt(tt2)[tt3] 
and 1 — a(l — 0:2) [1 — 0:3] giving rise to eight families of eigenf unctions. In order to get the 
ground state the parameters a, ^2, 03 should be chosen in such a way as to minimize the 
eigenvalue and then ( |0|) corresponds to the ground state. It can be denoted (Q;,a2,«3)- 
The other values of u, z/2, 1/3, if inserted in (|6.2|), describe (provided that the corresponding 
wavefunctions are normalizable) the ground states of the remaining seven other families 
of eigenstates: 

(a, ^2, 1 - as) , (a, 1 - ^2, 03) , (a, 1 - "2, 1 - "3) , 

(1 - a, a2, as) , (1 - a, 1 - ^2, 03) , (1 - a, "2, 1 - 03) , 

(1 - a, 1 - 0:2, 1 - 0:3) • 

So taking different values of u, z/2, z/3 in (6.2) one can exhaust all the types of ground states 
corresponding to the different families mentioned above.0 

Using the same approach as in Sections 2-4, we make a gauge rotation of ( |6.1| ) with the 

(s) —1 (s) 

ground-state eigenfunction as gauge factor, /i^cd = ~2\E'o '^bcd^o- ^ straightforward 
calculation leads to the operator (we omit an additive constant) 

AT 



h 



BCD 



= E + ^EM(^(a:. - x,))(a, - d,) + cot(l(x, + x,))(9, + d,)] 

i=l i<j 

N N ^, 

+ Z/2^COt(Xi)(9i + — ^COt(^)9i. 
1=1 ^ i=\ ^ 



(6.3) 



It is worth mentioning that if the operator ( |6.3| ) is written in the coordinates e*^, it appears 
in its rational form |]15| (cf. the rational form of the Ajv-i Sutherland Hamiltonian in 0). 

According to the above discussion about eight different families of eigenfunctions, it is 
sufficient to study the operator (6.3) for generic u, z/2, 1^3 and consider only eigenfunctions 
which are symmetric under refiections. In particular, this implies that eigenfunctions of 
(6.3) have the form 

py,V2,y3 ^ _P(coSXi, COSX2, . . . , COSXat) , (6.4) 

where F is permutationally symmetric. 

Next we encode the symmetry properties of the problem studied by introducing the 
permutation and refiection-invariant periodic variables (Tfc(cosa;), which satisfy 



N 



^ak{cos{xi))t = exp 

A:=0 



cos(xi))t" 



(6.5) 



^^See also the discussion in footnote 7 and in Section 5.1 

^•^In what follows, for the sake of simplicity, we continue to call ( |6.2|) the ground state, of course, 
keeping in mind above discussion. 
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(cf.(4.8), (5.11)). Let us emphasize that these variables (6.5) are characterized by the 
same period as the original Hamiltonian (6.1). In these variables the Hamiltonian h^BCj^ 
becomes 

N f) f) ^ f) 

"S5t.= E-4«5-a- + i:B.g^ (6.6) 



where 



Aij = iV(Ti_i(3-j_i -^[{i- l)cFi^lCFj+i + l)(Ti_/_iO-j+i_i 



l>0 

-(i - 2 - l)cri^2-i^j+i - {l + j + l)(3-i_i_iaj+i+i 

_^,(iV-j + i)(iV-j + 2)a,_2 , (6.7) 

and (jfc = 0, when k < ot k > N. The method of the calculation of the coefficients 
Aij,Bj is similar to that presented in App.C. This expression can be called the algebraic 
form of the BCn, Bn,Cn and Sutherland Hamiltonians. Similarly to what happened 
in all previously discussed bosonic cases the coefficients Aij, Bj are polynomials of second 
and first degree in ak, respectively. Hence, the Hamiltonian h^BCD t)e written in 
terms of generators of the Borel subalgebra of gl{N + 1) realized as first order differential 
operators (see Appendix A.l) as in Sections 2 and 5.1. The result is given in eq. (B.7) 
(see Appendix B). 

Then in accordance with the general definition given in [Q] we conclude that all 
BCN,B]y, Cn.Dn Sutherland models ( |6.1D are exactly solvable. The existence of the 
representation ( |6.6|) proves that there are infinitely-many eigenfunctions of the operator 
(5.3) having the form of polynomials in the variables ak- It also implies that totally 
(anti) symmetric- (anti) symmetric eigenfunctions with respect to permutations and refiec- 
tions of the BCjsi, Bj^,Cn and Sutherland models ( |6Tl|) have a factorizable form being 
the product of the ground-state eigenfunction (6.2) multiplied by a polynomial in the vari- 
ables These polynomials are related to finite-dimensional irreducible representations 
of the Lie algebra gl{N + 1) in the realization (A.l). They can be called the BCn Jack- 
Sutherland polynomials. 



6.2 The supersymmetric extensions 

The bosonic BCn,Bn, Cn and Sutherland models ( |6.1D have natural supersym- 
metric extensions. In this subsection we will construct these models and show that, as 
was the case for the other supersymmetric extensions discussed in Sections 3,4 and 5.2, 
these models are also exactly-solvable. Let us introduce the supercharges ( p.3| ) with the 
superpotential W given by 



1 1 ^ 

[log I sin(-(xi - Xj)) \ + log I sm{-{xi + Xj))\ + z/2 ^ log | sin(xi)| 



N 

W = 

i<j ^ ^ i=l 

N 



+ i^sl^logl sin(^)| . (6.^ 



i=l 
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(cf.(3.4),(5.8)). From the general formula ^^^bcd — 'we derive the supersym- 

metric Hamiltonians of the BCN,Bfyf, and D]^ Sutherland models: 

TV ^ 



'^sBCD 



1 d d 



+ ^— 1 



1 |., + 2.,-l+..A|]Uc 



(6.9) 



where C = ~^Y.f=i[i^{N — i) + U2 + As should be familiar by now, the next step 

is to introduce the gauge-rotated Hamiltonian h^^BCD = ~'^^g^'^'^sBcd^o- After some 
calculations we get 



N 



d 



d 



d 



d 



i=l 

+ 

N 

1=1 



,1 



dxi dxj 2 sin — Xj)) dOi 86 j 



COt(-(Xj + Xjjj{- h 



+ 0,) 



d d_ 

2--' ■ ^"""-dxi ' dxj' 2s\n^{\{xi + Xj)yd9i ^ Wj 



cot(x 



d 



d 



* dxj sin^fx,) 06 j 



^ i=i 



cot( 



9 



2^(9xi 2sin2(f)a^, 



(6.10) 



It is worth mentioning that if the operator ( |6.1CI| ) is written in the coordinates it 
appears in its rational form |^. 

As in Sections 3, 4 and 5.2 it is convenient to use the superspace formalism and to 
introduce the superspace coordinates Xi, which satisfies: 



N 



Y.Xk{cos{(f)i))t'' = exp 



fe=0 



E 

.n=l 



\n+l 



n 



-Sn{cOs{(f)i))f' 



(6.11) 



Here Xfe(cos(0j)) = (3-fc(cosXj) +a;Cfc(cosXj, —6i sinxj), where ak and (k are the elementary 
symmetric polynomials defined in ( |3.9| ) but of new arguments: x — > cos x,9 ^ —9 sin x. 
In terms of these variables the Hamiltonian (6.10) becomes 

^ N Q Q . N 

I dada ^ Aij 



h 



is) 

sBCD 



dXi dXj 



r ^ 
J daJ2^. 



d 



where 



Aj = Nxi-iXj-i - E [(^ - i)xt-iXj+i + (^ + J - 1; 



Xi-l'lXj+l-l 



(6.12) 



l>0 



B 



-{i-2- l)xi-2-iXj+i - i^ + j + ^)Xi-i-iXj+i+i + i[Xi+i-iXj-i-i - Xi+i-iXj-i-i 

+ Xi+l+lXj-l-3 — Xi+l+lXj-l-3 + Xi+l+lXj-l-1 — Xi+l+lXj-l-l] , 



1 -= f (j-N -l)x,-i- (1^2 + f + l + i^i2N- J -l))jxj 



.,y^N-j + l){N-j + 2)x, 



i-2 



(6.13) 



^"^This expression should be compared to the rational form of the supersymmetric Sutherland model 
(4.3) and the BCn Sutherland model (6.3). 



20 



The method of the calculation of the coefficients Aij, Bj is similar to that presented in 
App.C. In components h^sBCD becomes 



h 



sBCD 



where 



N 



N 

E 



-iy{N~j + l){N-j + 2) 



+ y + l + K2iV-J-l))j 



5 p a ■ 



i>0 



- (i - 2 - l)ai-2-iCj+i - {1+ 3 + l)<3"j-/-iO+«+i 

— /[o-j-/-l(Ci+/-l + Ci+Z+l) + 0"i+/+l(Ci-«-3 — Ci-«-i) 



l>0 



- {i-2- l)Q-2~lcrj+i - (l+j + l)Ci-l-iCrj+i+i 

+ ^[^j-l-l{Ci+l-l + Ci+l+l) + 0"i+«+l(Cj-Z-3 — Cj-/-l)] 



, (6.14) 



- (i-2 - l)Ci-2-iCj+i -{1+ 3 + l)Ci-/-iCj+;+i 

+ ^[Ci+i-lCj-Z-l — Ci+«+lCj-/-3 + 2Ci+/+lCj-/-l 



(6.15) 



lYiom. the expression ( |6.14[ ) with the coefficients ( |6.15| ) it follows, that h^sBco can be 
written as a quadratic polynomial in the generators of a Borel subalgebra of gl{N + 1|A^) 
and, hence, the BCn Sutherland model is exactly solvable. 

The existence of the representation of ( |6.14|) in terms of the generators of a Borel sub- 
algebra of gl{N + l|iV) (see (B.8)) proves that there are infinitely-many eigenfunctions of 
( |6.14|) having the form of polynomials. These polynomials are related to finite-dimensional 
irreducible representations of the algebra gl{N+l\N) in the realization (A. 3). These poly- 
nomials can be called the supersymmetric BCn Sutherland polynomials. 

So the supersymmetric BCn, Bn,Cn and Sutherland models possess an algebraic 
form and also the Lie-algebraic form (B.8) being represented by second-order polynomials 
in the generators of the of the algebra gl{N + 1|A^) with certain coefficients. 
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7 Conclusion 



In this paper we have described the rational, algebraic and Lie-algebraic forms of the in- 
tegrable Aj^ — Bj^ — C^ ~ D]^ rational (Calogero) and trigonometric (Sutherland) Hamil- 
tonians in addition to the superspace expressions for the supersymmetric generalizations 
of the these models. We have shown that 

All Hamiltonians of the integrable Ajy — B^ — CN — DN rational (Calogero) and 
trigonometric (Sutherland) models possess the same hidden algebra gl{N + 
1) and can be represented by second-degree polynomials in the generators of 
a Borel subalgebra of the gl{N + l)-algebra. If the configuration space is 
parametrized by permutationally symmetric coordinates v, the Hamiltonians 
have a triangular form and preserve the flag of spaces of inhomogeneous poly- 
nomials 

Vn = spaniv'l'vl^^v^' . . . t;]^^ |0 < J] < n} , 

Consequently, each Hamiltonian possesses one or several infinite families of 
polynomial eigenf unctions. 

All Hamiltonians of the supersymmetric generalizations of the A^ — Bn — 
Cn — Dn rational (Calogero) and trigonometric (Sutherland) models possess 
the same hidden algebra gl{N + 1|A^) and can be represented by second-degree 
polynomials in the generators of a Borel subalgebra of the gl{N +l\N)-algebra. 
If the configuration space is parametrized by permutationally symmetric coor- 
dinates f , K the Hamiltonians have a triangular form and preserve the flag 
of spaces of inhomogeneous polynomials 

Vn = span{<it;^^..t;]^'^/€™i/€^^..K^^|0 < ;^r2i + ;^mi < ri,mi = 0, 1} , 

Consequently, each Hamiltonian has one or several infinite families of polyno- 
mial eigenf unctions. The integrability of the Ajy supersymmetric system was 
proven in j^; as for the Bn — Cn — Dn systems, it is an open question. 

As an interesting issue for future study we would like to mention the question about 
whether there exists a Lie algebraic description of the higher A^ — B^ — Cn — D^ rational 
(Calogero) and trigonometric (Sutherland) Hamiltonians. It is very probable that this will 
be the case and, perhaps, even the Lax operator can be represented in terms of gl{N)- 
generators. There are almost no doubts that it will be possible to extend the analysis of 
Sections 2-6 to the case of the exceptional simple Lie algebras. In particular, a special 
case of the one-parametric ^(2) rational and trigonometric Hamiltonians corresponding 
to three-body interactions only ||22|, ^ possesses a hidden (7/(3)-algebra and, hence, is 



expressible in terms of the gl{?>) generators |2^ |^. 

The question about the integrability of the above supersymmetric models leads natu- 
rally to the question about the possibility to develop the Hamiltonian reduction method 
for the supersymmetric cases and attempt to connect these models to 2d supersymmetric 



gauge theories along the hues of [|Tl[]. 



^^where = 0. 

^^Recently, a Lie algebraic form for the general G(2) rational and trigonometric models was found 
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An alternative description of the supersymmetric models discussed in this paper can 
be obtained by using the matrix representation of the 6'"'s in terms of the Pauli matrices, 

(T" (g) 1 ® 1 (g) ■ ■ ■ 
(T+ ® 1 ® 1 ® ■ ■ ■ 
cr^ ® cr" ® 1 ® ■ ■ • 

® a+ ® 1 ® • • • 

(7.1) 

In this formalism our supersymmetric models become matrix generalizations of the 'scalar' 
A]\r — Bj^ — C]\f — D]\r Calogero and Sutherland models. For the case our matrix models 
are particular cases of the general construction with internal degrees of freedom proposed 
in |251. 

The existence of the Lie algebraic form for the exactly-solvable — — Cn — 
Calogero and Sutherland models can be considered as a good starting point to investigate 
whether there exist other exactly-solvable problems and/or quasi-exactly-solvable gener- 
alizations (for definitions see, for instance, [^, |2^). The first step in applying the Lie 
algebraic formalism to the search for quasi-exactly-solvable problems was taken in [p8| . 

As a final comment we would like to mention that a great challenge in the subject is 
the search for solutions of the A^-body elliptic model (for a description see, for example, 
IIIOII). For the two-body case, the rational and algebraic forms have been known since 
Hermite's days. Recently, it was found that the 2-body case (the Lame operator) admits 
a Lie algebraic form with the hidden algebra gl{2) and is a quasi-exactly-solvable problem 
p9[] (see also [|1|). It gives a hint that the three-body (and more generally the A^-body) 
elliptic problem has to be quasi-exactly-solvable if any analytic solution exists whatsoever. 
However, all attempts to find a rational, algebraic or Lie algebraic form even for the three- 
body elliptic model have failed so far. 
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Appendices 



A Representations of the Lie algebra gl{N) and the 
Lie superalgebra gl{N\M). 

A.l gl{N) 

The generators of the Lie algebra gl{N) can be reahzed by first order differential operators. 
For our purposes we need one of the simplest realizations of gl{N) acting on the real 
(complex) space of dimension (A^ — 1). This is the space of minimal dimension where this 
algebra can act. The generators can be represented in the following form: 

d_ 



Eoi = Ji = — , 2 = 2, 3, . . . , iV , 



d 

^ij = Jlj = '^iJj = ' ^, J = 2, 3, . . . , iV , 

N Q 

Eqo = = n-J^'^kT;— , 

Eio = Jt = riJ\ z = 2,3,...,N, (A.l) 



where the parameter n G R (C). One of the generators, namely J° + J2p=2'Jp,py 



IS 



proportional to a constant and, if it is removed, we end up with the algebra sl{N). The 
generators J^j form the algebra of the vector fields of sl{N — 1), which is a subalgebra 
of gl{N). If n is a non-negative integer, the representation ( [A.l| ) becomes the finite- 



dimensional representation acting on the space of polynomials in {N — 1) variables of the 
following type 

K(t) = span{rf Tg'^^rr ■ ■ -t^^IO <Y.n, < n} . (A.2) 

This representation corresponds to a Young tableau with one row and n blocks and is 
irreducible. The positive-root generators define the highest-weight vector. If the Jj^'s 
are removed, the remaining generators form a Borel subalgebra. 



A.2 gl{N\M) 

Similarly to what was done for the algebra gl{N) one can construct a representation of 
the Lie superalgebra gl{N\M) in terms of first order differential operators over the direct 
sum of a even space and an odd space. Again the simplest realization of gl{N\M) act on 
the space C**^^^^'*^^ spanned by the the even variables (r*|2 = 1, 2, . . . , A^ — 1)) and the 
odd variables {k'^I'J = 1,2, ... , M) and is given by 







-^00 = 




—n + S— 


+ E 

7 




Eio = J, 




Eoi = 


Jr = 


d 






E^i Qyi 




E- — 




r- ' 
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d 

Eo-y = Q~ = ^— ^ , -£"70 = = k'^T 

E,, = T^, = (A.3) 

where the parameter n G R (C). Here the indices run over the following values, i = 
1 . . . N — 1 and 7 = 1 . . . M. The subset {Eij, E^i, Ei^, E^p] form the "vector field" 
representation of gl{N — 1|M). The algebra has a Z2 gradation and it is clear that 
the elements with an odd number of Greek indices are odd and the elements with an 
even number are even. Using the composite index notation / = (0, i, 7), the commutation 
relations of the superalgebra can be written (cf. the commutation relations in the defining 
(A^ + M) X [N + M) supermatrix representation, with {Ejj)mn = ^in^jm) 

[Eij, Ekl] = SjkEil - SjlEkj , (both even) 

[Eij, Ekl] = ^jrEil - S^Ekj , {Eu even, Ekl odd) . (A.4) 
{Eij, Ekl} = SjkEjl + 6ilEkj , (both odd) 

The Lie superalgebra gl{N\M) is not simple. There are two different cases: (i) A^ 7^ M 
and (ii) A^ = M. Let us consider the first case, when A^ 7^ M. Removing the unit element 
Yl!i=Q ^ Eii we are left with the superalgebra sl{N\M), which is simple. The "diagonal" 
generators for this case are given by 

-^00 — Ell , 1 = 

J En - Ei+ij+i , 1 < I < N - 1 , . 

E11 + Ei^ij^i , I = N (^-^^ 

Eii - Ei+ij+i , N + 1<I<N + M~1 

When N = M the situation is slightly more complicated. There are two one dimen- 
sional abelian ideals, which must be removed in order to make the algebra simple. 

The generators of gl{N\M) defined in ( [A.3| ) act in the space spanned by the monomials 



{ti^ ■ ■ ■T^NSi'^i^ ■ ■ ■ i^^m}- Here 6k equals either or 1. When n (in the expression for 
above) is an integer we have finite-dimensional highest weight representations whose 
highest weight vectors satisfy Y.k=i + S^li 5fc = These representations form a flag 
of spaces. 

B The Lie algebraic forms of the A^-i^ BCn, Bn, Cn 
and Hamiltonians 

In this Appendix we collect all Lie algebraic forms of the Hamiltonians for the (super- 
symmetric) Ajsi-i, BCn, Bm, Cn and Calogero and Sutherland models found in 
and in the present paper. 

• An^i Calogero model 



The Lie algebraic form of the many-body Calogero model is the following [12 




25 



2 51 I ~ ^N^^ ^-Jj-i,jJk-i,k - ~ ^ + 2i)Jj+i_ijJk-e-i,k 



2<k<j<N K 



N / 1 \ ^ 

- 2^5]A;Jfe,fc- - + z/ 5](iV-fc + 2)(iV-A; + l)Jfc_2,fc (B.l) 

k=2 ^ k=2 

where the following notations are introduced for the gl{N) generators (A.l): Jij = 

TO jQ ^ J- jo =n 

The supersymmetric Calogero model has the Lie algebraic form: 



.(rel) _ f ( iN-t+l){j-l) 
"■sCal ~ \ AT [^i-l,j^j-l,i-r'Ji^l,iJ^j-l,j-tJj~i,jli-i. 



i>max(l,j— i) 

+ Jj-l-l,jTi^i-i^i — Tj+i_ijTj_/_i^j] > + ^^l[Tj^2~l,jJi+l,i — Ti^ijJj-2-l,i 



1>1 

N 

+ Ti^i^iJj_2-i,j ~ Tj^2~i,iJi+i,j — '^Ti^ijTj^2-i,i] ~ 2^^;^^ j[Jjj + Tj-j] 

i=2 

1 ^ 

- (^ + ^^)E(^-J + 2)(iV-J + l)[^,-2,+T,_2,], (B.2) 
i=2 

where the following conventions are used for the gl{N\N — 1) generators (A. 3): 
Jij = J°. = and Tij = T'°0, if i and/or j is less than 2 or greater than N. 

An -I Sutherland model 

The Lie algebraic form of the many-body Sutherland model in terms of the gl{N) 
generators is the following [O 



^-1 i j(N - ?1 ^ 1 

j=l { 1=1 J 1=1 



+ 2 E { ^^^^ ^^ J,,Jk,k-j2U-k + 2l)J,+i,Jk-i,k].iB.3) 

l<k<j<N-l { 1=1 ) 

Here we have used the identifications Jij = Jfj, Jo,i = Jf = JN,i, Ji,i = 0. The 
supersymmetric many-body Sutherland model written in terms of the gl{N\N — 1) 
generators has the form 

_ _ pi^ ~ r J. . J. , J. .r. . M T T T- T- 1 

''sSuth ~ \ l'^hj'^],t ^ '^hi'^ 3,] ^ 

+ E U ~ ~ + Ji+l,iTj-lJ + Jj-ljTi^i^i — Ti^ijTj^i^i 

(>max(l,j— i) 

+ ^^l-[Tj^l^ljJi+l-l^i — Ti^i^ijjj^i^ii + Ti^i^i^iJj^i^ij — Tj^i^i^iJi^i^i 
1>1 
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N-l 



i=l 



(B.4) 



In ( [B.4| ) Jjj = = and Tj = T-° = are by definition zero when i and/or j is 
greater than — 1 or less than 1. 

BCn, Bn, Cn and Calogero models 

The bosonic BCn, B^, Cn and Calogero models have the following Lie algebraic 
form in terms of the gl{N + 1) generators 



TV N 

hBCD = 4 ^ ^(2/ + 1+J- ^) J,_z_i,, J,+z,. - 4^ E 

ij=l l>0 i=l 
N . ^ 



+ 2J2^[1 + U2 + 2u{N - j)]{N - J + 1) - -(iV + 3 - 2z) j J,_i, (B.5) 



and the supersymmetric extensions have the form 

N 

hlicD = 2EE[2(2/ + 1 + j — i)[Ji-i-i,jJj+i,i + Ji-i-ijTj^i^i + Ti_i_ijJj^i^i 

i,j=l l>0 

~ Ti-i-ijTjj^i^i\ + l[Tj_i^ijJi^i^i — Ti^ijJj_i_i^i — Ti^i^ijJj_i^i 



N 

+ Ti+ijTj_i_i^, + Ti+i^ijTj^^]] + ^ 2[1 + z/2 + 2z/(iV - j)]iN - j + 1) 

i=l 

N N ^ 

.[-{N + 3- 2z)T,_i,, - (AT + 3 - 2t)J,_,A - Au Y.^[J^,^ + T,,,] . (B.6) 

in terms of the gl{N + 1|A^) generators. In the above two formulas, as well as in 
the two formulas below, we have used that Jij = J?- = and Tij = T^j = when i 
and/or j is greater than or less than 1. 

BCn, Bn, Cn and Dn Sutherland models 

Finally, the Lie algebraic form of the bosonic BCn, Bn, Cn and Dn Sutherland 
models is 



''BCD 



N 

E 



NJi^ijJj^i^i - E [(^ ~ l)Ji-i,jJj+i,i + + j - ^)Ji-i-i,jJ; 



l>0 



— {i — 2 — l)Ji-2~i,jJj+i,i — + j + ^)Ji-i-i,jJj+i+i,i 

N 

+ E 

i=i 



|(j - A^ - 1) J,_i,, - + I + 1 + AT + z.(2A^ - J - 



- [u{N-j + l)iN-j + 2) + 2N]J,^2,j 



(B.7) 
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in terms of the gl{N + 1) generators, and the Lie algebraic form of their supersym- 
metric extensions is 

N r 

Kbcd = X] ] N{Ji-i,jJj-u - - l)[Ji-ijJj+i,i + Ti-ijJj+i^i 

i,j=l I l>0 

+ Ji-i,jTj+i^i - Ti_ijTj+i^i] + + l)[Ji-i-ijJj+i-i^i + Ji^i-ijTj+i_i^i 
+ Ti_i_ijJj+i_i^i - Ti_i_ijTj^i_i^i] - (i - 2 - l)[Ji_2-i,jJj+i,i + Ti^2-i,jJj+i,i 
+ Ji-2-i,jTj+i^i - Ti_2-i,jTj+i^i\ - {1 + j + l)[Ji-i-ijJj+i+i^i + Ti_i_ijJj+i+i^i 

+ Ji-l-\,jTj+l+l,i — ri_;_ijTj+;+i^i] + /[Jj_i_ij(Tj+;_i^j + Tj+^+i^j) 

- Ji+i+i,i{'^j-i'-3,j - Tj-i-i,j) - Ti+i+ij{Tj_i_i^i - Tj_i_3^i) 



N 



{Ti+i-i,j + Ti+i_ij)Tj_i_i^i^ > + 

) i=i L 

1^3 



|(i-7V-l)[J._i,. + T,_i,,.] 



[u{N -j + l){N-j + 2)]( J,_2,, + T,_2,,) - 2iV( J,_2,, - r,_2,, 



fB. 



in terms of the gl{N + 1\N) generators. 

In conclusion, let us emphasize the fact that all Lie algebraic forms of the supersym- 
metric models (B.2), (B.4), (B.6), (B.8) contain no fermionic (odd) generators. 

C Deriving the Lie algebraic forms: technical details 

In this appendix we will give some details of the derivation of the Lie algebraic forms of 
the Hamiltonians discussed in this paper. We will exemplify the method for the BCn 
Calogcro models. The other cases are tackled in a similar way. 
(A) Bosonic BCn Calogero model. 

Our final goal is to derive (5.6)-(5.7). Let us recall that the gauge-rotated BCn Hamil- 
tonian is given by 



N 



i=l i<j 



^^id,-dj) + -^{d. + d,] 



N 



N 



+iy^Y,-di-2ujY,Xidi. (C.l) 



After the change of variables Xi di the Hamiltonian becomes 



f^i dxi dxi dak d^i ^ dxf ddk i^ii^i ^ 9xi dak 



28 



N 



]^—\ '^i ~l~ •^j dxj Xi Xj dxi dxj dcj^ 

^ 1 da, d 



' ^^~{Xi dxi dak 
As an example, let us express 



(C 



N 



^ ^ L -J* ■ —I— nr ■ riT ■ riT ■ T ■ — T* ■ ii'T ■ riT ■ 



and 



^ daidak 



j^—i dxi 



(C 



(C 



in terms of the a variables. The approach we will follow is to represent the generat' 
functions of (C.3) and (C.4) 

fe=0 fc=0 

in terms of the generating function of the a variables 



AT 



\n+l 



A;=0 



n=l 



n 



-S2„(x)r] . 



(C 



Let us begin by calculating the generating function B{t) for (|C.3| ) 

1 



B{t) = E 



d d . 1 
+ 



d d 
+ 



G(t) 



2E E(-i 

i^j ^ m=l 



,m+l 



2m-l I ™2m-l 2m-l _ 2m-l-, 



Xi ~\~ Xj 



r G'(t) . (c 



Using the relation 



2n-l I 2n-l J2n-1 _ 2n-l n-1 

+ ^ — = 2 2^ 



•3y 2 I T 



2Z 2(n-l)-2Z 
■^3 



and making rather obvious mathematical transformations, one can show that 



d 



At"— - A(2N - 2)t'^ + AtN{N - 1)J G{t) . 



Substituting (^) in (C.8) we find that B, in (|a3D is 

Bt = A{N - k + 1){N - k)ak-i . 



(C 



(C 



(C 
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Next we proceed to the calculation of ( |C.4| ). The corresponding generating function 
can be written as 



N 



dG{p) dG{t) 



N oo 



4EE(-i) 



fe+£^2(fc+^)-2,fc / 



X4 



t''p'G{t)G{p) 



dxi dxi j^j^ ^^^^-^ 
Changing the summation variables to r = A; + £ and u = k — i one obtains that 



(C.IO) 



N 00 



^(i,p) = 4EE(-ir^2(.-i)(tp)^E 



i=l r=2 

The sum over u is easily determined 



p 



G{t)G{p). 



(C.ll) 



(1 



1 - 



V 



(C.12) 



Using this formula it is easy to transform ( |C.11| ) to the following expression for the 
generating function: 

it rd d 



A{t,p) 



1-^ 

p 



i__p_JGWG(p). 



Expanding the rhs of (C.13) we arrive at the final expression for (|C.4|) 

q>0 



(C.13) 



(C.14) 



The other terms in ( |C.2[ ) can be calculated using the same method which finally leads to 

the results presented in (|5l6| ), (|5.71 ). 

(B) Supersymmetric BGn Calogero model. 

For the supersymmetric models the calculations are similar to those carried out for the 
bosonic models. In terms of the superspace ^-coordinates (3.11) the gauge rotated Hamil- 
tonian (|5.10| ) is written as 



>'sBCD 



dada E 



5 



+ 



■{ 6(f)i{a) 6(f)i{a] 
1 . S 



+ 2u 
5 



da^2 

i<j 



S 



(l)i{a) - (f)j{a) 50i(a) 5(j)j{a)' 



(l)i{a) + (j)j{a) 6(j)i{a) S(f)j{a) 

f ^ 5 
2u / da 4>i{(y) 



N 

+ z/2 da^ 



6 



\ (f)i{a) 6(f)i{a) 



1=1 



(C.15) 



Now we make the change of variables 



Xi (see (5.11)) using the relation (|3.21 



Finally, an expression similar to (|C.2|) emerges except for the fact that the bosonic coor- 
dinates are replaced by supercoordinates with integration over the Grassmann variables: 



N 

h^sBCD ~ E / dadadpdjS 
i=i 



(50,(a) 5<P,ia) Sxi{P) 5xiiP) 



30 



N 

2u / dadpY^Y. 



i<j k=l 



1 ^ Sxk{f3) ^ SxkiP) . 
(j)i{a) - (j)j{a) 6(j)i{a) S(j)j{a)' 
1 



(f)i{a) + (j)j{a) S(j)i{a) 6(j)j{a)' 



N N 



U2 J dadp ^ ^ 



1 SxkW) 6 
6(f)i{a) SxkiP) 



N N 



+ 2uj J (ja (f)i{a 



i=l k=l 



.(C.16) 



For terms linear in derivatives the process of calculation is exactly the same as for the 
bosonic case. The final result can be obtained from the bosonic calculation simply by 
replacing the cr-coordinates by x-supercoordinates. However, for the term quadratic in 
derivatives: 



N 



i=l 



/ dadadpdp 



sxkmsxm s 6_ 

S<P,{a) 5cP,{a) 5xiiP) Sxi{P) 



dpdpA!, 



kl 



6 6 
SxiW) SxiiP) 



(C.17) 



a slight complication arises. Again our strategy is to make use of a generating function 
of the X coordinates 



N 



GiP,p) = T.Xk{Mf3))p' = exp 



k=0 



■n=l 



n 



(C.18) 



in order to represent the generating function for ^^'^^ = YliLi I (iada %^{a) • 



N 



(C.19) 



k,e=i 



In what follows we use the notation 0j := and 0j := 

Let us first note that 



N oo 



+in i2n—l i2m—lj.n m 



(C.20) 



i=l m,n=l 



(cf.(C.lO)). After some simple mathematical transformations similar to those leading to 
(C.ll) A^'^^tjp) can be rewritten as 



2Yi-intpy^T.[^2r-2i^) + S2r-2 



r=2 
oo 



P 



'giP,p)g{P,t) 



r=2 u "^V ^) ^P^ 

Using ( |C.12 ) together with the relation 



u \P 



-(r-2)+r--r - + (r - 2) - 
p \pJ \p 



(C.22) 
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(^]2T|) takes the form 



r=l 



2-J--^T.-^M<P)-^2r{<l>)W-f)0if3,PmM . (C.23) 

p (1 — r) ~; r 



r=l 



Evidently, the first term in (|C.23| ) is equal to 



At 8 d 



(C.24) 



(cf. (|C.13|) ). In order to calculate the second term in (|C.23|) we need to use two relations: 

CO + l 



E 

r=l 



[S2M~S2MW-P" 



_y^°o (-1)' 

g Z—ir — l r 



-[s2r{4>)-S2rWW-Pn _ 1 



(C.25) 



and 



Finally, the second term in (|C.23|) becomes 

t t + p 



p{l 



L\2 



g{(3,t)G{(3,p)-G{(3,p)G{(3,p) 



(C.27) 



We should emphasize that the bosonic part of this expression vanishes. Combining (C.23) 
and (C.26) we obtain the final expression for the generating function A^^\t^p) in terms 
of the generating function Q{P,p): 



it , d d 



t t+p 



G{P,t)GCp,p)-G{P,p)GCp,p) 



(C.28) 



(cf.(C.13)). Expanding in powers of t^p we arrive at the explicit expression for 

the coefficients in front of the second derivative terms 



= 2 E ^ 2(2g + 1 + - ^)Xl-q-lXk+q - l[Xk+qXi-q-l " Xk+qXl^q-l 
q>0 



+ Xe+q-lXk-q — Xi+q-lXk- 



(C.29) 



(cf.(C.14)) given in ( |5.12| ). Making similar calculations for the coefficients in front of the 
terms linear in derivatives we come to the conclusion that these coefficients are equal to 
those given by (5.7) with replacement of a's by x's (see (5.7) and (5.13)). 
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